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ABSTRACT 

We report on observation results of the prompt X- and 7-ray emission from 
GRB011211. This event was detected with the Gamma- Ray Burst Monitor and 
one of the Wide Field Cameras aboard the BeppoSAX satellite. The optical coun- 
terpart of the GRB was soon identified and its redshift determined {z = 2.140), 
while with the XMM-Newton satellite, the X-ray afterglow emission was de- 
tected. Evidence of soft X-ray emission hues was reported by Reeves et al. 
(2002), but not confirmed by other authors. In investigating the spectral evolu- 
tion of the prompt emission we find the possible evidence of a transient absorption 
feature at G.Qtg^s keV during the rise of the primary event. The significance of the 
feature is derived with non parametric tests and numerical simulations, finding 
a chance probability which ranges from 3 x 10~^ down to 4 x 10~^. The feature 
shows a Gaussian profile and an equivalent width of 1.2^06 keV. We discuss our 
results and their possible interpretation. 
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1. Introduction 

The nature of the progenitors of celestial Gamma-Ray Bursts (GRBs) is a still open 
issue. Collapse of massive fast rotating stars (hypernova model, e.g., Paczynski 1998) or 
delayed collapse of a rotationally stabilized neutron star (supranova model, Vietri & Stella 
1998) are among the favoured scenarios for the origin of these events. Both models predict 
that the pre-burst environment consists of a high density gas, due to strong winds from the 
massive progenitor in the case of a hypernova or to a substantial enrichment of heavy elements 
by a previous supernova explosion (SN) in the case of the supranova model (Lazzati ct al. 
1999; Bottcher 2000; Weth et al. 2000). The recent discovery, from the spectral evolution of 
the optical afterglow emission from GRB030329, of a direct connection between this burst and 
the Type Ic energetic supernova (hypernova) SN2003dh (e.g., Stanek et al. 2003, Hjorth et al. 
2003) seems to point toward the hypernova model, even if the strict simultaneity of the two 
events has not been fully proved. It is however not yet clear whether all GRBs are connected 
with energetic super novae, lacking in many cases evidence of supernova spectroscopic features 
and/or a 'supernova bump' in the optical emission of GRB afterglows. On the other hand, 
from the current estimates of the opening angles of GRB jets, the expected GRB rates, 
when compared with those of SNe Type Ib/c, suggest that less than 1% of such SNe might 
be accompanied by a GRB (Frail et al. 2001). The case of SN2002ap (e.g., Wang et al. 
2003) is an outstanding example of Type Ic hypernova which has no simultaneous detected 
GRB event associated with it. 

The presence of a post-supernova environment can be tested from the study of the burst 
X-ray spectrum, which should show a low-energy cutoff and/or absorption features due to 
elements in the circumburst material left by the previous supernova explosion. Either the 
progressive photoionization of the neutral gas by the GRB photons (e.g., Bottcher et al. 
1999; Lazzati & Perna 2002) or the electron density decrease or the electron temperature 
increase of an already almost photoionized medium (Lazzati et al. 2001) would make the 
features transient. 

In fact a transient absorption feature at 3.8 ± 0.3 keV has been detected in the prompt 
emission X-ray spectrum of GRB990705 (Amati et al. 2000). The feature is present in the 
first, rising part of the burst profile and disappears thereafter. Interpreted by Amati et al. 
(2000) as a cosmologically redshifted K edge due to neutral Fe around the GRB location, 
a GRB redshift of 0.86 ± 0.17 was derived and later confirmed by optical spectroscopy of 
the associated host galaxy (zopt = 0.84, Le Foc'h et al. 2002). With this assumption, the 
Iron relative abundance with respect to the solar one was derived, Fc/Fcq = 75 ± 19, which 
is typical of a supernova explosion environment. An alternative explanation was given by 
Lazzati et al. (2001), who assumed that the feature is an absorption line due to resonant 
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scattering of GRB photons on H-like Iron (transition ls-2p, Erest — 6.927 keV). Also in this 

case the derived redshift is consistent with that of the host galaxy and the line width is 
interpreted as due to the outflow velocity dispersion (up to 0.1c) of the material, which 
should have a Fe relative abundance of ~ 10 with respect to the solar one. In both scenarios, 
the observed feature points to the presence of an iron-rich environment left by a recent 
supernova explosion. 

GRB011211 was localized with the BeppoSAX WFCs (Gandolfi 2001) and promptly fol- 
lowed up with both the XMM-Newton satellite and ground based optical/radio telescopes. In 
the optical band, a counterpart was soon discovered (Grav ct al. 2001; Bloom & Bcrger 2001) 
and its redshift determined {z = 2.140 ± 0.001; Fruchter et al. 2001, Holland et al. 2002). 
In the X-ray band, Santos-Lleo et al. (2001) reported the discovery of afterglow emission 
with XMM-Newton, while Reeves et al. (2002, 2003), analyzing the XMM-Newton data, re- 
ported the evidence of five fading emission lines which were consistent with blue shifted Ka 
hues from Mg XI, Si XIV, S XVI, Ar XVIII and Ca XX outflowing at a velocity of ~ 0.1c. 
However the significance level of this detection was questioned. Rutledge and Sako (2003), 
examining the same data analyzed by Reeves et al., found that the claimed lines would not 
be discovered in a blind search, showing, with Monte Carlo simulations, that these features 
would be observed in 10 percent of featureless spectra with the same signal-to-noise ratio. 
Also Borozdin & Trudolyubov (2003) questioned the reality of these features. 

In the context of a systematic investigation (Frontera ct al. 2004) of the spectral evolu- 
tion of all GRBs jointly detected with the Wide Field Cameras (WFCs, Jagcr ct al. 1997) 
and Gamma Ray Burst Monitor (GRBM, Frontera et al. 1997, Costa et al. 1998) aboard 
the BeppoSAX satellite, we have analyzed the prompt emission of this burst. Surprisingly 
we find a marginal evidence of an absorption feature which is visible only during the rise 
of the event. If true, this is the second case of a transient feature after that observed from 
GRB990705. In this paper we report the results of our findings and their possible interpre- 
tation. 

2. BeppoSAX observations 

GRB011211 was detected with the BeppoSAX WFC No. 1 and GRBM (Frontera et al. 2001) 
on 2001 December 11 at 19:09:21 UT. Since a follow-up with XMM-Newton was promptly 
scheduled for this GRB, no observation with the BeppoSAX Narrow Field Instruments was 
performed. The GRB position was determined with an error radius of 2' (99% confidence 
level) and was centered at q;2ooo = 11^15'"16^ and (Jaooo = -21°55'44" (Gandolfi 2001). 
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Data available from the GRBM include two 1 s ratemeters in two energy channels (40- 
700 keV and >100 keV), 128 s count spectra (40-700 keV, 225 channels) and high time 
resolution data (down to 0.5 ms) in the 40-700 keV energy band. The WFCs (energy 
resolution ^ 20% at 6 keV) were operated in normal mode with 31 channels in 2-28 keV 
and 0.5 ms time resolution. The burst direction was offset by 9.1° with respect to the 
WFC axis. With this offset, the effective area exposed to the GRB was ~ 500 cm^ in the 
40-700 keV band and 65 cm^ in the 2-28 keV energy band. The background in the GRBM 
energy band was estimated by linear interpolation using the 250 s count rate data before and 
after the burst. The WFC spectra were extracted through the Iterative Removal Of Sources 
procedure (IROS ^, e.g. Jager et al. 1997 ) which implicitly subtracts the contribution of 
the background and of other point sources in the field of view. 

Among the gamma-ray bursts simultaneously detected by the BeppoSAX WFC and 
GRBM, GRB 011211 is the longest in both X- (2-28 keV) and 7-rays (40-700 keV). Figure 1 
shows its light curve in these energy bands. In X-rays (top panel of Fig. 1) it shows a long 
rise (~ 300 s) and a shorter decay (~ 120 s). Secondary peaks are also visible: three during 
the GRB rise and one during the burst decay. In gamma-rays (bottom panel of Fig. 1), 
the count statistics do not allow a precise determination of the onset time, which seems to 
occur later (~ 100 s) than the onset of the X-ray emission. However the gamma-ray peak 
is achieved about 20 s earlier than in X-rays and the GRB end appears to be simultaneous 
in X- and gamma-rays. The GRB entire duration in gamma-rays is ~ 270 s. 



3. Spectral analysis and results 

Following the investigation performed on a sample of BeppoSAX GRBs (Frontera et al. 2000), 
we are performing the systematic analysis of all GRB detected with the BeppoSAX WFCs 
and GRBM, in order to study their spectral evolution in the 2-700 keV energy band. Our 
method is that of subdividing the pulse profile in a certain number of slices of duration 
such to derive statistically significant count spectra for each of them. For GRB011211, we 
subdivided the fight curve in 3 intervals A, B, C (see Fig. 1), 100 s, 76 s, and 116 s duration, 
respectively. The derived 2-700 keV count rate spectra are shown in Fig. 2. The surprise 
was the shape of the spectrum in the time interval A, which showed at a glance, between 6 
and 8 keV, a structured depression, which was not apparent in the other two spectra, which 
exhibited the typical shape found in many other GRBs. Subdividing the 100 s interval in 
two sub-intervals, the count statistics were very low, but, in spite of that, the depression 
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marginally appeared in both spectra. Subdividing the interval B in two sub-intervals, no 
trace of the feature was found in either spectra. 

The feature in the interval A is also evident (see Fig. 3) when we perform the ratio of 
the GRB count spectra with that of the Crab Nebula, which by chance was observed at an 
angular offset similar to that (9.1°) of GRB011211. The Crab ratio criterion, commonly used 
to search for cyclotron resonance features in X-ray pulsars (e.g., Dal Fiumc ct al. 2000), 
has the advantage of being independent of the detector response function adopted and thus 
minimizing the calibration uncertainties. 

We investigated a possible instrumental origin of this depression with negative results. 
The high-voltage, which is monitored every second, does not show any glitch in the interval 
A. Given the accuracy of the readout (which dominates over statistical noise) this excludes 
gain changes higher than 0.01%. All the ratemeters show the GRB event, even the ones that 
measure the illegal events, so there is nothing that can be concluded from them. There are 
no dips or spikes of any sort, with a typical 3a upper limit of 10% per second of measurement 
(for a count rate of 700 counts s~^). All monitored quantites at the time of the depression do 
not show anything out of the ordinary. Finally, note that the same procedure was followed to 
extract both the spectrum A and the spectra B and C, in which the feature is not observed. 

As for the other events, the joint WFC and GRBM count rate spectra A, B and C 
were analyzed by means of the XSPEC software package (Arnaud 1996). In the fits, allowing 
the normahzation factor of the GRBM data with respect to the WFC data free to vary, we 
tested various continuum models: power-law (pl), broken power- law (bknpl) and smoothly 

broken power-law (or Band law, BL, Band et al. 1993). The PL model provided a good 
fit only with a GRBM/WFC normalization factor of 0.3 which is unusually low. Indeed, 
from extended flight calibrations with celestial sources (e.g., Crab Nebula) and using GRBs 
observed with both the WFC-GRBM and the CGRO/BATSE experiment (e.g., Fishman et 
al. 1994), this factor was found to range from 0.8 to 1.3. Thus we excluded the PL model. 
The BL gave a good fit, but the best fit parameters derived were not constrained by the data. 
The best data description and the best parameter constraints were obtained with a bknpl, 
with best fit parameters reported in Table 1. In the fits the GRBM/WFC normahzation 
factor was fixed to 0.8, but we verified that the fit quality and the BKNPL parameter values 
did not significantly change by using other values of this factor in the 0.8-1.3 range. We 
also took into account, in the fits, photoelectric absorption of the radiation from the GRB, 
finding that the equivalent hydrogen column density A^h was consistent with zero. Thus we 
assumed a Galactic absorption along the GRB direction (Ah = 4.2 x 10^° cm~^. Dickey & 
Lockman 1990). 

Even if, for the interval A, the bknpl model already provides an acceptable fit (x^/dof — 
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20.4/22 corresponding to a chance probability of 0.57, see Table 1), the distribution of the 
residuals from the model (see Fig. 4) still shows evidence of the structured depression. The 
feature is still apparent in the residuals from the best fits with the other models (PL and bl) 
and is also evident when only the WFC data arc fit with a PL model. The used models arc 
those generally used to describe the prompt continuum emission from GRBs. However we 
cannot exclude that a properly tuned two-component continuum model could recover the 
feature, even if, on the basis of the results obtained thus far in the X-ray band (< 10 keV) 
with BeppoSAX and with HETE-2 (e.g., Barraud et al. 2003), this possibility appears much 
less probable. 

Thus wc decided to evaluate from a statistical point of view the significance of the 
feature using the residuals from the bknpl best fit model. Given that the x^-test is a global 
test, which does not take into account the relative position of the residuals, we performed 
a null-hypothesis test to check the independence of the residuals around the depression (see 
Fig. 4). For that we used the run-test (e.g., Bendat and Piersol 1971), a non-parametric test, 
which is valid independently of the statistics assumed for the data. Taking the residuals to 
the BKNPL as data points, we investigated the probability that groups of n consecutive 
data points are mutually independent. Conservatively the n value was chosen to be the 
number of energy bins corrersponding to 4 times the energy resolution of the instrument in 
correspondence of the depression region. Thus we considered values of n ranging from 11 
to 14 bins. For a chosen value of n, we scanned the entire count spectrum, by moving the 
first bin of the group. We found that the probability that n consecutive residuals from the 
BKNPL model are mutually independent, ranges from about 0.16 to 0.47 when the energy 
bins are far from the depression, and from 2.7 x 10~^ down to 7.2 x 10~^ (equivalent, in 
a Gaussian approximation, to 2.8-3.1 o"), depending on n, when the energy bins cover the 
depression region. The probability values obtained in the depression region for the different 
values of n are reported in Table 2. 

Stimulated and encouraged by the run-test results, we replaced the photoelectrically 
absorbed BKNPL model with the same law modified by a negative narrow Gaussian profile, 
according to the following equation: 



where a — cf{E) is the photo-electric cross-section of a gas with cosmic abundance (Morrison 
& McCammon 1983), and E^^ is the break energy, A is the normalization factor, EW is the 
equivalent width of the line, and 
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is the assumed Gaussian profile, with centroid energy Ei and standard deviation ai. A similar 
model is, e.g., used to describe cyclotron absorption features from X-ray pulsars (e.g., Soong 
et al. 1990). The fit results with this new model, instead of the bknpl, are reported in 
Table 1. In the fits ai is set to 0.1 kcV. Leaving it free to vary, the best fit is still obtained 
for (T/ ^ 0.1 keV with an upper value of 1.2 keV at 90% confidence level. In Fig. 5 wc show 
the A spectrum with the best fit curve and the residuals of the data from the model. As can 
be seen, the depression, has now disappeared (x^/dof = 12.1/20). 

Being aware of the warnings given by Protassov et al. (2002) about the use of the 
F-tcst for the inclusion of an additional term to a fitting function (Bcvington 1969) when 
this function is an absorption or emission feature, in order to evaluate the probability that 
the parent fitting function of the spectrum A is actually given by the Eq. 1, we followed 
the following approach. Assuming as featureless spectrum the bknpl model that best fits 
the WFC plus GRBM data (see Table 1), by means of xspec we generated 10 000 WFC 
fake spectra with the same number of energy bins of the measured A spectrum. These 
spectra take into account the time duration of the spectrum A, the response function of 
the instrument and the Poisson noise associated to the measured spectrum. We fit each 
of these spectra with a bknpl model, obtaining a set of 10 000 distributions of residuals 
from the model assumed. For each of these residual distributions, we performed a blind run- 
test, determining the spectra which have at least one group of n consecutive residuals from 
the BKNPL whose probabihty of being mutually independent is less than the probabilities 
reported in Table 2. The fraction of the spectra that satisfy this condition are reported 
in Table 3, for values of n in the range from 11 to 14, as above. 

Given that the run-test takes into account the relative position of the residuals but 
not their amplitude, we also fit these spectra with the Eq. 1 model, and, mindful of the 
Protassov et al. (2002) warnings, on these spectra we performed not the F-test for testing 
the addition of a further component to a fitting function, but for testing the discrepancy 

between an assumed fit function and the parent function. This test (see Bevington 1969, p. 
195) makes use of the probability distribution of the ratio F12 = xli/xlz ('^^ -^21 = l/-^i2) 
where xt-^ = X^/^i is obtained from the fit of the function 1 (in our case the BKNPL model) 
to the data, and — is obtained from the fit of the function 2 (in our case, Eq. 1). 

In the spectral fits with the Eq. 1, the ai was again frozen to 0.1 keV, and, in order to 
perform a blind search of a negative Gaussian feature, for each spectrum we performed 6 fits, 
each with a different starting value of the feature centroid. In this way we determined the 
fraction of spectra which satisfy the condition F12 > F^, where is the value measured 
of F12, derivable from the values reported in Table 1. We also checked the condition 
F21 < -^21) finding the same result. The values of ft versus n are reported in Table 3. As 
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can be seen, the fraction of simulated spectra which satisfy both tests ranges from 3 x 10~^ 
down to 4 X 10~^. These fractions give our best estimate of the statistical significance of the 
negative Gaussian feature observed. 

In the above estimates we did not take into account the number of trials whereby such 
a feature could be found as a result of a random fluctuation. Thus, assuming that the 
range of values found for ft as chance probability range for a random fluctuation, using the 
binomial distribution, we evaluated the probability P(2) that, observing n GRBs, two of 
them (GRB990705 and GRB011211) show during their leading edge an absorption feature. 
To determine n, from the entire sample of GRBs detected with both WFC and GRBM we 
selected those which show during their leading edge 2-10 keV count spectra with statistical 
quality at least similar to that of the spectra in the time intervals in which the absorption 
features are apparent (interval A for GRB011211, interval B for GRB990705). We found 
n = 18 and thus 2.4 x 10"^ < P(2) < 1.3 x 10~^, depending on the chosen value of the 
chance probability p, 4 x 10~^ or 3 x 10~^, respectively. The low values of P{2) gives us an 
additional confidence that the feature observed is not a random fluctuation, at least on the 
basis of the GRB sample size obtained with BeppoSAX . 

Assuming the best fit centroid energy of the negative Gaussian and ai — 0.1 keV, we 
also evaluated the 2a upper limit to the feature EW in the time intervals B and C. These 
upper limits and the corresponding are also reported in Table 1. 

We have also tested a K edge model (edge in XSPEC) in the place of the negative 
Gaussian. In this case the fit improvement is lower, likely due to the symmetrical shape of 
the line which does not favor this model. 

On the basis of the best fit results, we have derived the energetics associated with 
the event. The fiuence of the GRB is 5'2_io = 1-1 x 10~^ erg cm""^ erg cm~^ in the 2-10 
keV band, and 5*40-700 = 5.1 x 10~^ erg cm~^ erg cm~^ in the 40-700 keV band, while 
the total 2-700 keV fiuence is 5*2-700 = 7.5 x 10~^ erg cm~^. For comparison, in the time 
interval A (100 s duration), the 2-700 keV fiuence is ~ 2.6 x 10~^ erg cm~^. The peak 
flux averaged on a 8 s bin is -Ff-io ~ ^-^^ x 10^^ erg cm~^ s~^ in the 2-10 keV band, and 
-^40-700 ~ ^-0 ^ 10~^ erg cm~^ s~^ in the 40-700 keV interval. /^From the optical redshift, 
assuming a standard cosmology with Hq ~ 65 km s~^ Mpc~^, fl\ = 0.7 and Qm = 0.3, and 
isotropic emission, the total fiuence corresponds to a released energy of ~ 3.55 x 10^^ erg, 
with 34% of it released during the first GRB interval. The released energy in the time 
interval A above 7 keV (feature centroid energy) is 1.2 x 10^^ erg, which corresponds to 98% 
of the energy released by the burst in this interval. 
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4. Discussion 

There is a long history of X-ray hnes from GRBs. In the eighties, evidence of hard 
X-ray (>10 keV) absorption features was reported by Mazets et al. (1981), Hueter (1984), 
and Murakami et aL (1988). The most significant result was obtained by Murakami et 
al. (1988) with the Gamma Burst Detector aboard the Ginga satellite, with two transient 
absorption hnes observed during the latest part (5 s duration) of the GRB880205 leading 
edge. The line energies were ~ 19 and ~ 38 keV and their chance probability was evaluated 
to be 6 X 10~^ (Fenimore et al. 1988). Another significant hard X-ray line detection, which 
was visible, like the previous one, during the leading edge of the burst, was later detected 
with Ginga from GRB890929 and reported by Yoshida et al. (1991). 

All these features at those times had a natural explanation in terms of cyclotron fea- 
tures from Galactic neutron stars (GNS). Once this association with GNS resulted to be 
untenable owing to the isotropic distribution of GRBs discovered with BATSE, the general 
belief was that all these feautures likely were not real. Possible justifications included in- 
accurate knowledge of the true response functions of the instruments used, partly due to 
the indetermination of the GRB direction. Another tacit justification was the difficulty of 
making these lines at cosmological distances. No detection of hard X-ray lines was reported 
by the BATSE team. 

With the advent of the afterglow era, evidence of X-ray lines has been found both 
in the GRB afterglow and in the prompt emission. Concerning the afterglow, emission 
lines with energy centroids below 10 keV, have been reported from several GRBs (see, e.g., 
review by Frontera 2003), most of which pointing to Iron fluorescence or recombination 
lines from circumburst regions. Concerning the GRB prompt emission, as discussed in 
Section 1, a transient absorption feature was detected during the rise time of GRB990705 
(Amati et al. 2000), In that case, the identification of the absorption feature, either with 
a K-edge of neutral Iron (Amati et al. 2000) or with a resonant scattering line of nearly 
ionized Iron (Lazzati et al. 2001), led to a redshift coinciding with the host galaxy redshift 
independently measured in the optical band (Le Floc'h et al. 2002). This fact strenghtened 
the result. 

Now we have the possible evidence of another similar feature from GRBO 11211. Also 
in this case, the feature is visible only during the rise of the primary event. However, in 
the case of GRBO 112 11 the interpretation is not so straightforward. First of all, the feature, 
due to its symmetrical shape, is more consistent with an absorption line than a K edge. 
Also the assumption of a cyclotron resonance feature, which was the assumption for the 
origin of the Ginga lines, appears untenable. If an oriented magnetic field were present 
in the emitting region, the magnetic field strength can be derived from the condition that 
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the magnetic energy density integrated over the emission region of the prompt emission 
times a radiative efficiency factor (hkely not smaller than ~ 0.1) is at most equal to the 
estimated energy released (see Section 3). The result is that B = 10^-10^ Gauss. This 
appears inconsistent with the magnetic field strength (~ 10^^ Gauss) obtained assuming 
that the line is a cyclotron resonance feature {Erest = 21.7^J g keV after correcting it for 
the optically-determined redshift). Neither this large magnetic field seems to be compatible 
with the value derived from the analysis of the afterglow data (Jakobsson et al. 2003) using 
the standard homogeneous fireball model (Sari et al. 1998). 

The rest-frame energy of the absorption feature is also inconsistent with any plausible 
resonant scattering line unless we invoke a blue shift of the line. Indeed the most abundant 
elements of typical interstellar environments or those typical of a supernova explosion have 
resonant scattering energies far below the rest-frame energy of the discovered feature (see, 
e.g., Kato 1976). This implies an absorbing medium with very high outflowing velocities 
Vout- A lower limit to this velocity is obtained assuming a resonant scattering of the GRB 
photons off H-like Ni XXVIII (rest-frame energy of 8.092 keV). To derive the outflow speed, 
we consider that photons seen in the onflow comoving frame aX v = vq will be absorbed, 
where z/q is the resonance frequency of the considered line. The following condition needs 
therefore to be satisfied: 

{1 + Z)6 Z/abs = t'O (3) 

where z/abs is the observed absorption frequency and 5 = [7(1 — f3cos6)]~^ is the Doppler 
factor. In our case, since the photon and the fluid direction coincide, we have cos 6' = — 1 
and the equation can be rewritten, after simple algebra, 

(1 + ^)^ (4) 



where Hvq is the rest frame energy of the candidate feature. In this case we get (5 — Vout/c — 
0.75, with a corresponding Lorentz factor 7 = 1.5. Note that in this trans-relativistic regime 
nor the /? ~ 1 neither the 7 ~ 1 approximations can be adopted. Eq. 4 simplifies though to 
the usual formula z) z^abs = 271/0 in the ultra-relativistic regime. 

Assuming, as proposed by Lazzati et al. (2001) for GRB990705, that the absorption 
feature is due to resonant scattering of GRB photons off a mixture of H-like Fe and Co 
{Erest = 6.927 and 7.518 keV, respectively; oscillator strength of 0.416, Kato 1976), the best 
fit (x^/dof = 13.1/20) is obtained assuming a super-thermaP Gaussian distribution of the 



thermal distribution of outflow velocity would give a too narrow resonant feature, since fcTpe < 10 keV 
in order to have a sizable fraction of non fully ionized iron ions. 
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product (5^ where 7 is the bulk Lorentz factor of the expanding medium, with up^/ {(3^) — 
7.5% and neghgible Fe K-edge of the absorbing medium. The best fit parameters of the 
model (mean value of /S'-f, and column density of the Fe absorbing material) are reported in 
Table 1. From the best fit value of f3'-f, we derive, by using again Eq. 4, a value of j3 = 0.80 
and 7 = 1.65. The deconvolved spectrum of the interval A with superposed the best fit 
model is shown in Fig. 6. The low depth of the Fe K edge could be the result of the high 
ionization level and/or the edge broadening due to Compton scattering (Ross et al. 1999). 
The disappearance of the feature after ~ 100 s is likely the result of the inhibition of the 
electron recombination caused by the increase of the electron temperature as a consequence 
of Inverse Compton heating (Lazzati et al. 2001). 

The high outflowing velocity is a problematic issue. The acceleration cannot be due to 
the interaction with the burst photons themselves, since the feature centroid does not evolve 
when the interval A is analyzed with higher temporal resolution. A possible explanation is 
that the absorption feature is due to outflowing material in the line of sight toward the jet 
cone (see also Konigl 2004). 

The fit with the physical model allows us to estimate the distance of the absorber and 
the amount of absorbing material. Prom Eq. 12 of Lazzati et al. (2001) we derive: 

i?Abs - 5 X 10^^ cm (5) 

while the fit yields A^fc+Co — 10^° cm~^. Assuming isotropy, the mass of the absorber results 
in Mpe+co = O.I5M0, the total mass being M ~ {60/Ape)MQ where Ape is the iron (cobalt) 
richness in solar units. This estimate implies a moderate iron and cobalt enriched material 
in order to keep the total mass to a reasonable value. Alternatively, one can envisage 
a clumped material, with a covering factor much lower than unity, with a dense clump 
lying randomly along the line of sight. In addition, it is possible to derive the total energy 
involved in the outflow, from the best-fit average speed. We obtain E'kin = (7 — l)Mc^ ~ 
5 X 10^^ /[A-Fer]{Q/ATr)] erg where rj is the covering factor of the clumps and Q is the solid 
angle within which the material is accelerated to the high speed obtained from the fit. 
Reasonable values are Ap^ ~ 10, rj ~ 0.1 and f2/47r ~ 0.1, resulting in a large amount of 
kinetic energy E'kin ~ 5 x 10^^ erg, a value larger than the typical GRB energy of ~ 10^^ erg 
(Bloom et al. 2003). 

Given that any absorption feature in the GRB prompt emission, specially during its 
rise, implies that the absorbing material is external to the fireball, there is only one model, 
among the proposed progenitors, which can accommodate this observation. In the Supranova 
scenario (Vietri & Stella 1998) a supernova explodes prior to the GRB, ejecting high speed 
material around a newly formed rotating supra-massive neutron star. The GRB is produced 
when the spin of the neutron star is reduced and the compact object collapses to a black 
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hole. In the time span between the supernova and GRB explosions, the spindown of the 
neutron star releases E^b ~ 10^^ erg, of which a significant fraction could be in the form of a 
pulsar-type wind. As discussed in Konigl & Granot (2002), this energy might be deposited 
in a pulsar-wind bubble and used to accelerate the supernova shell. They in fact argued that 
the swept-up remnant is likely to expand anisotropically, with the polar regions of the shell 
(which can be identified with the gas responsible for the high-energy absorption detected 
in GRB 011211) attaining speeds ~ c.^ On the other hand, explaining the presence of a 
large amount of high velocity absorbing material in the immediate vicinity of the burster is 
extremely problematic for any model invoking a simultaneous explosion of the burst with 
the progenitor star (e.g. Woosley 1993), since the stellar wind creates an iron poor (solar at 
most) low-density environment at the required distance. 

Only more sensitive X-ray observations than those performed with the BeppoSAX 
WFGs can enlarge the GRB sample, improve the statistical quahty of the time resolved 
spectra, and thus test the results obtained with BeppoSAX on the X-ray transient absorp- 
tion lines. These observations are expected to be only partly covered by the upcoming 
satellite Swift (e.g., Krimm 2004 given that these lines are observed only during the leading 
edge of the bursts, while the Swift Burst Alert Telescope has an energy passband from 15 
to 150 keV. However Swift can solve the issue of the hard X-ray lines observed with Ginga. 
A sensitive mission which would be capable to observe the GRB since their onset and to 
test the presence of transient X-ray lines with centroid energies below 10 keV could be the 
proposed Lobster-ISS for the International Space Station, an all-sky monitor with focusing 
optics from 0.1 to 3.5 keV (Eraser et al. 2002) paired with a GRB monitor for higher photon 
energies, from 2 keV to several hundreds of keV (Prontera et al. 2003b; Amati et al. 2004). 

We wish to thank Nicola Masetti and Eliana Palazzi for useful comments. This research 
is supported by the Italian Space Agency (ASI) and Ministry of Instruction, University and 
Research (COFIN funds 2001). BeppoSAX was a joint program of ASI and the Netherlands 
Agency for Aerospace Programs. 
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Table 1. Spectral parameters of the best fit models 



Parameters A B C 



BKNPL 





1 




1 cin+0.21 

i.<JU_g 45 


r2 




1.69l°«, 


2-26l^;2i 


Eh (keV) 


12.41^^0° 




13lf 


xVdof 


20.4/22 


22.1/27 


16.4/27 


BKNPL modified by a negative Gaussian 


(Eq. 1) 






1 9+0.4 
-^•^-0.8 


0.951?;^? 


-L-OU.Q 47 


r2 


-■-•O-O.l 


1.69^°:°? 




Eb (keV) 




6-9^2:7 


istf 


El (keV) 


D-y-0.5 


[6.9] 


[6.9] 


EW (keV) 


1 9+0.5 
-■-•^-0.6 


< 0.7 {2a) 


< 0.9 (2a 


xVdof 


12.1/20 


21.7/26 


16.3/26 


BKNPL plus a H-like Fe 


resonance feature (see text). 












r2 








Eb (keV) 


7 2"'"^'^ 






A^FeXXVI+Co XXVII (cm~^) 


9.4l|^ x' 1019 








-1 09+0.16 
-'-•'-'^-0.16 






xVdof 


13.1/20 
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Table 2. Probability Pr, obtained with the run-test, that n consecutive data points 
around the region covered by the depression feature are mutually independent. 



n 


Pr 




11 


2.1 X 10" 


-3 


12 


1.2 X 10" 


-3 


13 


7.2 X 10- 


-4 


14 


2.7 X 10" 


-3 



Table 3. Fraction of simulated featureless spectra fr which provide a number of runs with 
probability P < Pr (see Table 2) and fraction of spectra ft which satisfy both the above 
run test condition and the condition on the F-test (see text). 



11 0.0262 0.0030 

12 0.0103 0.0009 

13 0.0044 0.0004 

14 0.0244 0.0021 
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Fig. 1.— Light curve of GRB011211 in two energy bands, 2-28 keV (WFC), 40-700 kcV 
(GRBM), after background subtraction. The zero abscissa corresponds to 2001 Dec. 11, 
19:04:31 UT. The time shces on which the spectral analysis was performed are indicated by 
vertical dashed lines. 
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Fig. 2. — Count spectra of the burst in the time intervals A, B, and C. 
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Fig. 3. — Ratio of the WFC spectra in the time interval A, B and C with the Crab Nebula 
spectrum measured by the WFC with the source in same offset as that of GRB011211. 
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Fig. 4. — Top: Count spectrum of GRB011211 in the interval A. The step-hke curve repre- 
sents the best fit with a bknpl. Bottom: Distribution of the residuals to the model. 
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Fig. 5. — Top: Photon spectrum of GRB011211 in the time interval A. The step-hke curve 
represents the best fit with the model given by eq. 1. Bottom: Distribution of the residuals 
to the model. 
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Fig. 6. — Top: photon spectrum of GRB011211 in the time interval A. The step-Uke curve 
represents the best fit with the model used in Lazzati et al. (2001) (see text). 



